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We study the particle number production and its time variation using non-equilibrium quantum field theory. We study the 
model proposed by Hotta et.al. [T| for particle number production with a heavy neutral scalar and a light complex scalar . The 
interaction Lagrangian contains CP violating phase and particle number violating interaction among the scalars. The particle 
number violating mass term is also introduced, which splits a complex scalar into two real scalars with small non-degenerate 
mass. Therefore, the term generates particle and anti-particle mixing. We study the long time behavior of the particle number 
production rate. 
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1 Introduction action. 

Study of the mechanism of the particle number pro- 
duction is a very important issue in baryogenesis and 
leptogensis. In experiments, there are many phenom- 
ena which violate the particle number. Such phenom- 
ena includes BB mixing and neutrino flavor oscillation. 
While they can be treated with time evolution of pure 
state, the baryogenesis and leptogenesis occur in the 
environment where the statistical treatment is suitable. 
This is because they occur when the definite particle 
number of the universe is unknown. In this study, We 
use the non-equilibrium field theory with the density 
matrix and study the time evolution of the particle 
number. This paper is organized as follows. In sec- 
tion II, we propose a particle number violating model 
which consists of a heavy neutral scalar and one com- 
plex scalar. In the next section, the current associated 
with the particle number is written in terms of a Green 
function of non-equilibrium field theory. In section III, 
the particle number production rate is computed and 
its property is discussed. The final section is devoted 
to conclusion. C = 
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N 



N 



B 2 (/> 2 +AN(f> 2 + h.c, 



(1) 



where TV is a real scalar and (f> is a complex scalar. 
There are two types of the interaction. One is parti- 
cle number conserving interaction which coefficient is 
given by and the other is particle number violat- 
ing interaction with the coefficient A. There are also 
two types of mass term. One of them with the coef- 
ficient B 2 violates the particle number and the other 
one is a particle number conserving one given by the 
mass term m 2 ^ '</>. One may take a phase convention 
that B 2 is real and A is complex. We denote the 
phase A as <f>A = arg.A and it is a source of CP vi- 
olation. The mass term B 2 breaks U(l) symmetry and 
it splits one complex scalar fields into the two mass 
eigenstates of real scalars. Introducing two real scalars 
as 4> = -^(01 + i4>2~)i the Lagrangian is rewritten as, 



^{d^Nd^N -M 2 N N 2 ) 
1^ 



2 ^ 



Lagrangian for the scalar model and Particle 
Number Production 



with 



two matrix and is given by, 



ij 

f B 2 . Aij is a two by 
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In the previous paper pQ , we proposed the following A = 
model for particle number production with the inter- 
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3 Computing the current for U(l) charge; par- 
ticle number 

The particle number associated with the complex field 
4> is a U(l) current. 



(4) 



jfj, = ity'd^ - 

= <hdn<t>\ -dft<h<l>i- 



We compute the divergence of the U(l) current with 
some initial condition specified with the density ma- 
trix. 



(i„(X))=Tr(j M (X)p(0)) 



(5) 



where p(0) represents an initial quantum statistical 
state; 



P(0) 



exp[-/3(g - fiL)] 
Tr[exp[-/3(ff - fiL)}} 



(6) 



L is U(l) charge given by L = J d 3 xjo and /3 is y with 
the temperature T. n is chemical potential. We as- 
sume that at t = all the interaction terms including 
U(l) breaking mass term are zero and t > 0, suddenly 
they are switched on. Therefore at t=0, Hamiltonian 
-ffo and particle number L commute as; 



[H ,L] = 0, 



(7) 



and at later time t > 0, the U(l) breaking terms are 
switched on and the particle number is not conserved. 



[H,L]?0. 



(8) 



In t = 0, Eq.([7|). the density matrix can be written by 
the following product, 



P(0) 



exp(—/iL) exp(— /3Hq) 
Tr[exp(-AiL) exp(-/3iJ )] ' 



(9) 



With the property Eq.0, one can write the density 
matrix even for non-zero chemical potential case. Us- 
ing the density matrix, the initial particle number at 
t = is given as, 



(i(0)) 



TrL(0)p(0) 
Vd 3 k 



sinh /z/3 



(27r) 3 cosh f3ujk — cosh p/x ' 



(10) 



When chemical potential is zero, the initial particle 
number is zero. In the following, we use the density 
matrix of /J, = and focus on the particle number pro- 
duction with the interaction. 



4 Green function and U(l) current 

U(l) current defined in Eq.([5]) can be written in terms 
of the Green function of non-equilibrium field theory; 

G\ 2 2 (x,y) = <Mx)p(o)Mv)> 

= Tr(0 2 (y)<MxV(O)). (11) 

The lower indices distinguish the species of the light 
scalar fields; cf>i(i = 1,2). The upper indices distin- 
guish whether the operator is on the time ordered path 
1 or on anti-time ordered path 2 in closed time path 
formulation [5]. Using the definition, the divergence 
of the averaged current with the density matrix is the 
production rate per unit time and unit volume. This 
can be related to the Green function as follows, 



d 



{J^X)) = (a x -D y )G\ 2 2 (x,y) 



x=y=X 



(12) 



If the space translation invariance holds, the current 
depends on only time. If this is the case, the diver- 
gence of the current is equal to time derivative of the 
particle number density, 

^<^)> = ^<J°PO>. (13) 

With the help of 2 PI (Two Particle Irreducible) for- 
mulation of the non-equilibrium quantum field theory 
[6] , we derive Schwinger Dyson equation for the Green 
functions in pQ. The Schwinger Dyson equation can 
be solved perturbatively and the divergence of U(l) 
current obtained in one-loop level. We find, 

d -(J^PO) ~ -4|A|sin< 



d 3 p 



d 3 k 
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2w p (27r) 3 J 2w fe (27r) 3 2w N 
{(n k + l)n N (n p + 1) - n k (n N 

{I(u 2 p + L02k - un,X°) 

-I(u lp + wife - u> N , X")}. 



l)n p } 



where we have ignored the terms which is proportional 
to B 2 and bJik = \/mf + k 2 (i 



(14) 
onal 

1,2). The time depen- 



dent function / is given as, 

cos flX° 



I(Q,X L 



(15) 



where njy ,nk(n p ) denote the thermal equilibrium dis- 
tribution function for N and <f> given as, 

1 

n N 



exp(/3wAr(fc + p)) — 1 ' 
1 



n k = 



exp(/3w fc ) - 1 



(16) 
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In the divergence of the current, first we study the fac- 
tor related to distributions, 



Substituting the time dependent factor, the divergence 
of the current is, 



(n k + l)n N (n p + 1) - n k (n N + l)n p . 



(17) 



The first term implies the decay N to two light scalars 
while the second term implies inverse decay. If the en- 
ergy conservation in the following sense, holds, 

uj k + uj p — uj N {p + k), (18) 



{J"{X)) ~ -4|A|sin^ A A 



d 3 p 



d 3 k 



2w p (27r) 3 J 2w fe (27r) 3 2ia N 
{(n k + l)n N (n p + 1) - n k (n N + l)n p ] 

d2 td2 



(22) 



then, 

{n N + l)n k n p = n N (n k + l)(n p + 1). (19) 

The particle number production is cancelled between 
the decay process and inverse decay process. Therefore 
the net particle number can not be produced if the en- 
ergy conservation of Eq. (|18p is satisfied. Next we study 
the coupling constants in Eq. (fl4|) . The production rate 
is proportional to CP violating phase (ft a ■ In addition 
to CP violation, both types of the interactions, one 
is CP conserving and particle number conserving one 
A$ and the other is CP violating and particle number 
violating one A are required. Finally we study time 
dependence. We are interested in the time length for 
which coherence of the two amplitudes in Eq.ffTJJ cor- 
responding to the N — s> <fii(k)<fii(p) (i — 1,2) is not 
lost yet. When the coherence remains, one can use the 
approximation, 



I(UJ 2P + L0 2 k - UJN,X°) - I(u>l p + UJ lk - U) N , X°) 

(20) 



amn X° . B 2 B 2 . „. 



' 2tu k 2co p ' 



where Oq = — w k — uj p . We have used the approxi- 
mation, 



0J2k = Jk 2 +m 2 +B 2 



UJ k 



B 2 
2uj k ' 
B 2 



u Jlk = ^k 2 + m 2 -B 2 = w fc -— . 

We consider when X° is large and is proportional to 
■ga-. In the limit of small B 2 with fixed B 2 X°, the time 
dependent factor becomes, 

I(UJ 2P + 0J2k - LO Nl X°) - I(u>l p + UJ lk - 0J N ,X°) 

(21) 



-Wsin { (^ + |-)X0 } 



The presence of the delta function implies the energy 
conservation ujn = u k + lj p holds. Therefore the decay 
contribution and inverse decay contribution are can- 
celled each other. Below, we consider only the decay 
contribution. 
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-4tt\A\ sin <PaA,i 
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( e Puj k _ i)( e ^ fc _ i)( e /3wjv - i) 

r>2 r>2 

xsin{(^ + ^)X°} 
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\A\A^T 2 sm<j) A 
v mam (u) 



sinh ~ Jo 



dV- 



U 2 -V 2 



cosh -j — cosh y 



where 



U = /3(u k + lj p ), V = /3(u k -u p ), 



and 



Vmax — \ (1 



4m3, 



in 



(23) 



(24) 



(25) 



We introduce the following dimensionless quantities, 
f3 = Pm 4> ,X = m 4 ,X ,B = —,m N = —. (26) 



To study the long time behaviour of the production 
rate X° ~ 4j , we define the rescaled time t, 



t = 



B 2 X° 



B 2 X { 



(27) 
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Figure 1: Time dependence of the production rate of the particle number density. The dashed line corresponds 
to $ = jq, the thick solid line corresponnds to $ = ^ , and the thin solid line corresponds to $ = ^. 



Using Eq. (|27|) , one may write the time dependent part 
as, 



dU 



$ 2 J/3m N sinh ^ 



2-irpUt 



rV ma AU) sm . 

o cosh 2 — cosh 



(28) 



In Fig.l, we show the time dependent factor Eq. (|28|) 
of the production rate as a function of the dimen- 
sionless rescaled time t for three inverse temperature 
$ = jq, 5g>arid -L. We choose rriN — 20. We can see 
the the production rate oscillates and the amplitude 
decreases. As temperature grows, the production rate 
becomes larger. 

5 Conclusion 

• We compute the time variation of the particle 
number with a scalar model. 

• In the model, a neutral scalar and one complex 
scalar are included. U(l) charge related to the 
complex scalar is the particle number. 

• The particle number and CP symmetry are vio- 
lated due to the mass term and the interaction. 

• Due to U(l) soft-breaking term, one complex 
scalar is not a mass eigenstate and mass eigen- 
states are two real scalars with non-degenerate 
mass. 

• We have computed the time dependence of the 
production rate of the particle number density. 

• The rate oscillates and is damping in their am- 
plitude. The oscillation period is larger for the 
smaller non-degeneracy. 



• Because of the damping, the particle number pro- 
duced during the first half cycle of the oscillation, 
may remain after the integration of the rate with 
respect to time. 

• We have not included the effect of the expand- 
ing universe, which effect leads to the decay and 
inverse decay contribution may not cancel. 

• As an extension of our work, one can consider 
the case for < L(0) >^ and how the ini- 
tial particle number is washed out. Also we 
may apply the method to flavored leptogenesis; 
L e (t),L^{t),L T (t). 
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